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Spin-orbit coupling (SOC), the intrinsic interaction between a particle spin and its motion, is 
responsible for various important phenomena, ranging from atomic fine structure to topological 
condensed matter physics. The recent experimental breakthrough on the realization of SOC for 
ultra-cold atoms provides a completely new platform for exploring spin-orbit coupled superfluid 
physics. However, the SOC strength in the experiment, determined by the applied laser wavelengths, 
is not tunable. In this Letter, we propose a scheme for tuning the SOC strength through a fast 
and coherent modulation of the laser intensities. We show that the many-body interaction between 
atoms, together with the tunable SOC, can drive a quantum phase transition (QPT) from spin- 
balanced to spin-polarized ground states in a harmonic trapped Bose-Einstein condensate (BEC). 
This transition realizes the long-sought QPT in the quantum Dicke model, and may have important 
applications in quantum optics and quantum information. We characterize the QPT using the 
periods of collective oscillations (center of mass motion and scissors mode) of the BEC, which show 
pronounced peaks and damping around the quantum critical point. 
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SOC plays a major role in many important condensed 
matter phenomena and applications, including spin and 
anomalous Hall effects [l| , topological insulators Q , spin- 
tronics f^, spin quantum computation, etc. In the past 
several decades, there has been tremendous efforts for de- 
veloping new materials with strong SOC and new meth- 
ods for tuning SOC with high accuracy for spin-based 
device applications [4,, i^J. However, the SOC strength in 
typical solid state materials (e.g., ~ 10"* m/s in semicon- 
ductors) is generally much smaller than the Fermi veloc- 
ity of electrons (~ 10^ m/s), and its tunability is limited 
and inaccurate. 

On the other hand, the recent experimental break- 
through on the realization of SOC for ultra-cold atoms 
provides a completely new platform for exploring SOC 
physics in both BEC |7Hl2j and degenerate Fermi gases 
P^16j . In a degenerate Fermi gas, such SOC strength 
can be at the same order as (or even larger than) the 
Fermi velocity of atoms. Because of the strong SOC, 
spins are not conserved during their motion and new ex- 
otic superfluids may emerge. For instance, new ground 
state phases (e.g., stripes, phase separation, etc.) may 
be observed in spin-orbit coupled BEC |8l4l2j and new 
topological excitations (e.g., Weyl and Majorana 13 1 
fermions) may appear in spin-orbit coupled Fermi gases. 
The observation and application of these exciting phe- 
nomena require tunable SOC for cold atoms. Unfortu- 
nately, the strength of the SOC in the experiment Q and 
other theoretical proposals [l7l - l20j is not tunable because 
the SOC strength is determined by the directions and 
wavelengths, not the intensities, of the applied lasers. 



In this Letter, we propose a scheme for generating tun- 
able SOC for cold atoms through a fast and coherent 



modulation of the Raman laser intensities [2l]|, which 
can be easily implemented in experiments. Such tun- 
able SOC for cold atoms provides a powerful tool for 
exploring new exotic Bosc and Fermi superfluid phenom- 
ena. Here we focus on a quantum phase transition (QPT) 
[2^ in a harmonic trapped BEC induced by the many- 
body interaction between atoms and the tunable SOC 
strength. With the increasing SOC strength, there is a 
sharp transition for the ground state of the BEC from a 
spin balanced (i.e., equally mixed) phase to a spin fully 
polarized phase beyond a critical SOC strength (i.e., the 
quantum critical point). By mapping the spin-orbit cou- 
pled interacting BEC to the well-known quantum Dicke 
model (23i . i24j . we obtain analytic expressions for the 
quantum critical point and the corresponding scaling be- 
haviors for the QPT, which agree well with the numeri- 
cal results obtained from the mean-field Gross-Pitaevskii 
(G-P) equation for the BEC. 

The realization of QPT in the Dicke model using the 
spin-orbit coupled BEC opens the door for many signif- 
icant applications in quantum optics, quantum informa- 
tion, and nuclear physics ^25.-27j . Previously the Dicke 
model has been studied in several experimental systems 
[H, [2§] , especially atoms confined in an optical cavity. 
However the coupling between atoms and optical cavity 
fields is very weak, and the experimental observation of 
the QPT in the Dicke model only occurred recently using 
the momentum eigenstates for a BEC confined in a cav- 
ity {sell . Compared with the cavity scheme, the spin-orbit 
coupled BEC utilizes the many-body interaction between 
atoms and has the advantage of essentially no dissipation, 
fully tunable parameters, very strong coupling, and the 
use of atom internal states, thus provides an excellent 
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Figure 1: (Color Online) An illustration of the experimen- 
tal scheme for realizing SOC for cold atoms (a) Laser 
setup, (b) Atom laser coupling, (c) A typical density distri- 
bution of one spin component of the BEC in the ground state, 
(d) Corresponding momentum distribution. The vertical and 
horizontal dotted lines are kx = and fcy = respectively. 



platform for exploring Dicke model related applications. 

Finally, the QPT is characterized using collective os- 
cillations of the BEC, such as the center of mass (COM) 
motion and the scissors mode, where the oscillation pe- 
riods show pronounced peaks at the quantum critical 
point. Furthermore, the oscillations of the BEC have 
regular periodic patterns in both spin balanced and po- 
larized phases, but show strong damping in the transition 
region. 

System and Hamiltonian: The harmonic trapped BEC 
in consideration is similar as that in the recent bench- 
mark experiment [6]. For simplicity, we consider a two- 
dimensional (2D) BEC in the xy plane with a strong 
confinement (with a trapping frequency uj^) along the z 
direction. Such 2D setup does not affect the essential 
physics because the z direction is not coupled with the 
SOC. Two hyperfine ground states jt) = \F — l,mp — 
— 1) and 14,) = \F = l,mF = 0) of *^Rb atoms define the 
spins of atoms, which are coupled by two Raman lasers 
(with Rabi frequencies fli and incident at a 7r/4 angle 
from the x axis, as illustrated in Fig. la and lb. 

The dynamics of the BEC are governed by the nonlin- 
ear G-P equation 

ihd^/dt = (pV2to + V{r) + Hs + Hi) $, (1) 



under the dressed state basis |t) = exp (iki • r) |t), |i) — 
exp (ik2 • r) jj,), where ki and k2 are the wavevectors of 
the lasers. $ = ($-[-, $4,)^ is the wavefunction on the 
dressed state basis and satisfies the normalization condi- 
tion / dxdy{\<^>^\^ + |$;|2) = 1. 
potential (r) = iTOw^(77^x'' 
ping frequency in the y direction, and 77 — ujx/ujy 
ratio of the trapping frequencies. Hs = jPxO'z + hfla^/^ 
is the coupling term induced by the two Raman lasers 
with az and (Jx as the Pauli matrices. The SOC strength 
7 — hkL/m, kL = |ki — k2|/2 — v^tt/A, and A is the 
wavelength of the Raman lasers. The Raman coupling 
constant = riif72/2A with A as the detuning from the 



The harmonic trapping 
"^), where ujy is the trap- 
is the 



excited state. The mean field nonlinear interaction term 

where the inter- and intra-spin interaction constants 
AttH^N^cq + C2)/'inaz and 



.9tt 



9n 



9U 



AttH Nco/muz, cq and C2 describe the corresponding s- 
wave scattering lengths [3l|, N is the atom number, and 
ttz = y^2Trh/mujz- 

Because the SOC strength 7 is determined by the laser 
wavevector kL, the SOC energy can be comparable to or 
even larger than other energy scales (e.g., the Raman 
coupling n) in the BEC. In a Fermi gas, 7 can be larger 
than the Fermi velocity of atoms. Unfortunately, due to 
the same reason, 7 cannot be easily adjusted in exper- 
iments, which significantly restricts the applications of 
the SOC in cold atoms. 

Tunable SOC for cold atoms- We propose a scheme for 
tuning the SOC strength 7 through a fast and coherent 
modulation of the Raman coupling $7 = fio + f7cos(a;t) 
that can be easily realized in experiments by varying the 
Raman laser intensities [sl]. Here the modulation fre- 
quency Lu is chosen to be much larger than other energy 
scales in Eq. ([IJ. In this case, the Hamiltonian in Eq. 
([T]) can be transformed to a time-independent one using 
a unitary transformation ip — exp[zr2 sin(wt)cra;/(2a;)]$. 
After a straightforward calculation with the elimination 
of the fast time- varying part in the Hamiltonian [s^, [sj] , 
the nonlinear G-P equation ([T]) becomes 

ihdip/dt = [p^/2m + V{r) + Hs+ Hi] ^, (2) 

where the Raman coupling becomes Hg — ^effPxCz + 
hfloax/^, with the effective SOC strength 



7e// = jJoi^l/uj)- 



(3) 



Here Jq is the zero order Bessel function. Clearly, 7e// 
can be tuned from the maximum 7 without the modula- 
tion to zero with a strong modulation. The mean field 
interaction term Hi = adV'tP + IV'J.P)^' + /^PV', a — g^^, 
P — ~ fftt) /2 s-iid P is a 2 X 2 matrix whose elements 
aregivenbyPn = |VtR|-^o(S) + iJo(f )] + l^iRi- 
i^o(^)], r22 = WVi - \M'-^)] + Wil + ^o(S) + 

lM^)h andPi2 = P^i = -i[l-Jo(f )][V'|V'i-^t^:]. 

We choose the physical parameters to be similar as 
those in the experiment 0: {ujy,ujz) — 2Tr x (40,400) 
Bz^r] = 1, X = 804.1 nm, cq 100. 8605, C2 = -0.46aB 
[351 with the Bohr radius as, N ^Ix 10*, uj = 2t: x 4.5 
kHz. For the numerical simulation, we need a dimension- 
less G-P equation that is obtained by choosing the units 
of the energy, length and time as hojy, ^Jh/ (rriLOy) — 
1.7 /im, and l/iOy = 4 ms, respectively. The dimen- 
sionless parameters in the G-P equation become 7 ~ 
^h/{mujy)kL = 9.37, a = 27V ^2 Trmujz/fi{co + C2) = 495 
and /3 = - N ^/2Trrnuj^C2 = 1.14. 

Quantum phase transition: The tunable SOC, in com- 
bination with the many-body interaction between atoms. 
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Figure 2: (Color Online) Quantum phase transition with tun- 
able 7e//. fio = 16, 7 = 9.37 is the bare SOC strength with- 
out modulation, (a) Plot of the spin polarization |{(Tz)| and 
{(Ta;) in the ground state. The blue and red lines are from the 
prediction of the Dicke Hamiltonian. The circles and squares 
are from the numerical simulation of the G-P equation ((2|. 
(b) Plot of the COM motion period T. The shift of the har- 
monic trap D = 1, To = 2-n /ujy. 



can drive a quantum phase transition between different 
quantum ground states in a harmonic trapped BEC. Here 
the ground state of the BEC is obtained numerically 
through an imaginary time evolution of the G-P equation 
. A typical density profile of the ground state is shown 
in Fig. Ic, which has a Thomas-Fermi shape, similar as 
that in a regular BEC. However, the momentum distri- 
bution of the BEC has a peak around the single parti- 
cle potential minimum located at {Kx, Ky) — (— iiTmin, 0) 
(see Fig. Id), where X„ii„ = \Jllff - ^o/^Te// and the 
degeneracy between ztii'njin is spontaneously broken. 



To characterize the ground state of the spin-orbit cou- 
pled BEC, we calculate the spin polarization \{crz)\ = 



and {a^) — 2Re/ dr^fjl^fj^. Here 



we choose the absolute value of (cr^) because the two de- 
generate ground states at iXmin have opposite (ctz) due 
to the spin-momentum locking term PxUz and they are 
spontaneously chosen in experiments. In Fig. 2a, we 
plot \{crz)\ and (ax) with respect to Jeff- For a small 
Jeff, the spin up and down atoms have an equal popu- 
lation, thus {<7z) = 0, {(Jx) = —1, i.e., the spin balanced 
phase. Beyond a critical point 7ejy, the spin population 
imbalance increases dramatically and reaches the spin 
polarized phase (tr^) = 1 {{(Jx) = 0) within a small range 
of "feff. The spin balanced and spin polarized phases 
at small and large jef / can be understood from the sin- 
gle particle Hamiltonian H^, where the Raman coupling 
flo<7x/'2 and SOC jeffPx^z dominate at the small and 
large jef f respectively. More numerical results show that 
the critical transition point occurs at "f^ff 



The QPT can be understood by mapping the spin-orbit 
coupled BEC to a quantum Dicke model. For an interact- 
ing BEC in a harmonic trap with a large atom number N 
(so that the mean field theory works), all atoms are forced 
to occupy the same many-body ground state. Therefore 
the energy variation for the change of the spin (e.g., spin 
flip) of one atom need be determined by the coupling 
between the atom spin and the many-body ground state 
mode. This is very different from an non-interacting BEC 
where atoms do not affect each other, but the same as 
that for many atoms interacting with a single photon 
mode in a cavity |23j . Treating the interacting many- 



body ground state as a single mode composed of different 
harmonic trap modes, we can map the Hamiltonian for 
the spin-orbit coupled BEC to 



H - hcUxNa^a + nn^Sx + 2^^lI^L^^{a^ - a){S+ - 5_), 

v2 

(4) 

which is similar to the Dicke model for two-level atoms 
coupled with a cavity field [2^. S^^y^z are the large spins 



for all atoms, 



Sy -j— %S z, S — 



iSz, a^a is 



a harmonic trap mode, a = y^mu!x/'2h{x + ipx/muJx). 
The critical point for the QPT can be derived from the 
standard mean-field approximation [23l, yielding the re- 
lation 7gjy. = ^fio/2, which is exactly the same as 
that from numerically simulating the G-P equation 
Just beyond the critical point leffi the Dicke model 
predicts that the scaling of the order parameters is 
\{az)\ = 2|^,|/iV = ^1 - {llffheff)\ {<y.) = 2SxlN = 

-{llffhefff iOTJeff > llff, and 1(0-^)1 = 0, {(Tx) ^ -1 

for Jeff < leff- Such scaling behaviors are confirmed in 
our numerical simulation of the G-P equation (see Fig. 
2a). The perfect match between numerical results from 
the G-P equation and the predictions of the Dicke Hamil- 
tonian shows the validity of the mapping to the Dicke 
model. 

We emphasize that the many-body interaction between 
atoms plays a critical role in the QPT by forcing all atoms 
in a single mode. Without interaction, numerical simula- 
tion of the G-P equation shows {az} = in certain region 
of 7e/ / > 7g^y , which disagrees with the prediction of the 
Dicke model. This disagreement indicates that atoms in 
a non-interacting BEC do not response to the change of 
Jeff collectively, although non- interacting and interact- 
ing BECs share the same transition for the energy spec- 
trum at 7e// 1 which changes from one single minimum at 
Kx = to two minima at ibfCmin. While for interacting 
BECs with large atom numbers iV = 4 x 10'* and 10^, we 
obtain exactly the same results as that in Fig. 2a, which 
further confirm the validity of our mapping to the Dicke 
model in the large N limit. 

Collective dynamics in BEC: the signature of QPT: 
It is well-known that various physical quantities may 
change dramatically around the quantum critical point 
(i.e., critical phenomena), which provides additional ex- 
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Figure 3: (Color Online) Plot of the COM motion period T 
versus 7e// and f2o- The insets show the corresponding l^x (t)). 
The green line is the theoretical prediction j^ff ~ \/0q/2 
from the Dicke model. The circles are the numerical results 
from the G-P equation Tq ^ ^n/ujy, D = 1. 



perimental signatures of the QPT. We focus on two types 
of collective dynamics of the ground state of the BEC: 
the COM motion and the scissors mode induced by a 
sudden shift or rotation of the harmonic trapping poten- 
tial, respectively. In a regular BEC without SOC, the 
COM motion is a standard method to calibrate the har- 
monic trapping frequency because the oscillation period 
depends only on the trapping frequency [s^ and is not 
affected by other parameters such as nonlinearity, shift 
direction and distance, etc. 

We numerically integrate the G-P equation ^ and cal- 
culate the COM (r(i)) = J dxdy{\^p-^{t)\^ + \^i^{t)\^)r (t). 
The COM motion strongly depends on the direction of 
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Figure 4: Plot of the scissors mode oscillation period T with 
respect to Jeff for three different Qq- 21 = 2-k jujy^J'rf- -f 1 
is the oscillation period without SOC, r\ = \/5, 6 = 4°. The 
inset is a typical density distribution of the ground state of 
the BEC. The dotted lines show the corresponding quantum 
critical points predicted by the Dick model. 



the shift T) of the harmonic trap. When D is along the y 
direction, the period of the COM motion along the y di- 
rection is Tq = 27: /ujy and not affected by 7e//, while the 
COM motion in the x direction disappears (i.e., {x) = 0). 
Here the COM period T is obtained through the Fourier 
analysis of (r [t)). The physics is very different when D 
is along the x direction, where {y (t)) = as expected, 
but {x (t)) depends strongly on Jeff, as shown in Fig. 
2b. In Fig. 3, we also plot T as a function of jeff and 
rio. Without SOC (7e// = 0), T = To, the period for a 
regular BEC, as expected. T increases with jeff in the 
spin balanced phase, but decreases when spin starts to be 
polarized, leading to a sharp peak at the quantum critical 
point '^eff- "^^^ oscillation of {x {t)) in the spin balanced 
phase is completely dissipationless, while a strong damp- 
ing occurs in a small range of jef f beyond j^^j: (see the 
inset in Fig. 3). Far beyond 7^^-^ , the oscillation becomes 
regular again with the period T — Tq because the ground 
state has only one component in this region. The peak 
and the damping of the oscillation around j^^j: provide 
clear experimental signatures for the QPT. Moreover, T 
also depends on the magnitude D of the shift near the 
critical point J^jf- the larger D, the smaller T. 

Another collective dynamics, the scissors mode [37|, 
shows a similar feature as the COM motion. The scis- 
sors mode can be excited by a sudden rotation of the 
asymmetric trapping potential (i.e., 77 ^ 1) by an angle 
9, which induces an oscillation of the quantity (xy) — 
J dxdy{\^^{t)\'^ + \ipi{t)\'^)xy. Without SOC, the period 

of the scissors mode is Ti = 2n/ ^J^x + o^" 
served in experiments [11]. In Fig. 4, We plot the oscil- 
lation period T with respect to Jeff for three different 
51o- We have confirmed that the same QPT occurs for 
\{<Jz) \ and {ux) of the ground state in this asymmetric po- 
tential with the quantum critical point j^j^ = -^/ilo/S, 
as predicted by the Dicke model. Similar as the COM 
motion, we observe the peak and damping of the oscilla- 
tion around llff Far beyond 7e//, the oscillation period 
is Ti . Similar as the dependence of the COM motion on 
the shift distance D, the angle 9 also influences the pe- 
riod of the scissors mode near I'^ff- the smaller 6*, the 
larger T. 

In summary, we show that the SOC strength in the re- 
cent breakthrough experiment for realizing SOC for cold 
atoms can be tuned through a fast and coherent modu- 
lation of the applied laser intensities. Such tunable SOC 
provides a powerful tool for exploring spin-orbit coupled 
superfluid physics in future experiments. By varying the 
SOC strength, the many-body interaction between atoms 
can drive a QPT from spin balanced to spin polarized 
ground states in a harmonic trapped BEC, which real- 
izes the long-sought QPT in the Dicke model and may 
have important applications in quantum information and 
quantum optics. 
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